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Unintegrated gluon distribution of a nucleus and photon-jet correlations 

Wolfgang Schafer a * 

a Institute of Nuclear Physics PAN, ul. Radzikowskiego 152, 31-342 Krakow, Poland 

An impact parameter dependent unintegrated gluon distribution is constructed as a solution of a nonlinear 
evolution equation with realistic Glauber-Gribov input. Photon-jet correlations in pA collisions in the proton 
fragmentation region are proposed as a direct probe of the nuclear unintegrated glue. 
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Dipole scattering amplitude vs. 
grated glue 



uninte- 



At high energies, where the relevant transverse 
momenta of partons fulfill Aqcd ^ P± \/s, 
we should take them explicitly into account. The 
most appropriate description is one in terms of 
unintegrated parton (dominantly gluon) distribu- 
tions. The latter are related to observables by the 
so-called fcj_-factorization. If we are in a regime 
where multiple scattering/saturation effects are 
large, such as in the case of a nuclear target, the 
linear fcj_ factorization breaks down and observ- 
ables are in general nonlinear functionals the un- 
integrated glue [l|2j . The nuclear unintegrated 
glue <p(b,x,p) is defined in terms of the forward 
scattering amplitude of a qq dipole r at impact 
parameter b, T(b, x, r): 

^-r(b,x,r)e-*P r = (l-w Q (b,x))6^(p) 



(27T)2 



(b,x,p) , 



(1) 



where wo(b, x) has the meaning of Bjorken's gap 
survival probability. To construct the nuclear un- 
integrated glue, one starts at moderately small 
x ~ xa ~ 0.01, where only the qq state is coher- 
ent over the whole nuclear size. Then the dipole 
amplitude for the nuclear target takes the simple 
Glauber-Gribov form: 



r(6, x A , r) = 1 - cxp[-a(x A , r)T A {b)/2] 



(2) 



where Ta(6) is the nuclear matter density. The 
nuclear unintegrated glue can then be expressed 
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as an expansion over multiple convolutions of the 
free-nucleon unintegrated glue f(x,p) [3]: 



<i>(b,x A ,p) = ^2w {b,x A )f {j) {x A ,p) 
Here 



(3) 



is the collective glue of j overlapping nucleons, 
and 



Wj(b,XA) 



v A (b,x A ) 



exp[-VA(b,x A )] 



is the probability for j nucleons to contribute. It 
depends on the effective opacity 



(4) 



v A {b. x A ) = 7;a s (q 2 ) cr (x A )T A (b) , 

where oq{x) — j d 2 pf{x,p) is a nonpertur- 
bative parameter, the cross section of a large 
color dipole. Interestingly, the expansion ((3|) 
also has a simple unitarity cut interpretation: 
namely <p(b, x, p) is proportional to the quasielas- 
tic quark-nucleus cross section, and the j-th term 
in eq. (|3]) is the contribution from j cut Pomerons 

01- 

2. Nonlinear evolution of the unintegrated 
glue 

At smaller values of x the multigluon-Fock 
states qqg, qqgig 2 ■ ■ ■ qqgi ■ ■ ■ g n with strongly or- 
dered light-cone momentum fractions z n <C 
z„_i ... « zi « 1 must be accounted for. 
While for the free nucleon target the effect of all 
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Figure 1. Unintegrated gluon distribution 
p 2 4>{vA,x,p) for different opacities va- Top: 
x = 0.01, lower: after evolution at x — 10~ 6 . 



multigluon Fock states is summed up by the linear 
BFKL-equation, no such simple procedure exists 
for the nuclear target. Simply iterating the first 
step of small-x evolution, one would arrive at the 
Balitsky-Kovchegov (BK) [5] equation. In terms 
of <f> it can be put in the form 



d<t>{vA-,x,p) 
d\og{l/x) 



= K-BFKL ® (fiiVA i X, p) 

+ Q[4>]{va,x,p), 



(5) 



with a linear, BFKL, piece, and a quadratic 
"gluon fusion" term Q[(f>]. The latter has a "real" 



as well as a "virtual" piece oc <t>{vA, x ,p)- The nu- 
merical treatment is much simplified by the obser- 
vation, that the real piece can in fact be written 
as a square: 

J d 2 qd 2 K(t)(vA,x,q) K(p + n,p + q) 

-K(p, K + p)- K(p, q + p) 4>{va,x, k) 



d 2 K 4>{va, x, k) 



p 



P + K 



P 



(P- 



2 



and can be viewed as a diffractive cut of a triple 
Pomcron contribution. At large p 2 , above the 
saturation scale, we obtain for the fusion term 



Q,[<j>]{v A ,x,p) oc — \ 



p 



d 2 K 



,2 



(VA,X,K 



\va,x 



P 2 ) I ~~Y ( d 2 qct)^A,x,q 2 ) 
Jp 2 K Jk 2 



Notice that it is pure higher twist, and involves 
only the "anticollinear" integration domain k 2 > 
p 2 . In particular it cannot be written in terms of 
the square of the integrated gluon distribution. 

In the practical solution of the evolution equa- 
tion eq.© one integrates over all transverse mo- 
menta, from soft to hard, and some regularization 
in the infrared domain is inevitable. For example, 
we freeze the running coupling as at small mo- 
menta. It is also necessary to introduce a finite 
gluon correlation radius (Xq 1 , to remove unphysi- 
cal long-range gluon exchange contributions. No- 
tice that this can be done without upsetting the 
gauge cancellations in the kernel, as it would be 
the case e.g. for a sharp momentum cutoff. 

Results from a numerical solution are shown 
in fig[TJ The boundary condition at xa = 0.01 
was constructed in terms of a free-nucleon glue 
fitted to HERA data. We used fi% = 0.5 GeV 2 . 
We plot p 2 (f){vA,x 1 p) ex dG/dp 2 , which is propo- 
tional to the phase space density of gluons. In- 
stead of impact parameter b we changed variables 
to the effective opacity va- Larger va corresponds 
to more central, and smaller va to more periph- 
eral collisions. 

The position of the maximum in fig[T]is a good 
definition of the saturation scale. We observe, 
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that it is a function of the opacity as well as of 
x. It increases for more central collisions (where 
absorption is stronger) as well as for smaller x. 
Notice that the evolution of the saturation scale 
with va is a result of the dynamics, and does not 
involve an assumption its impact-parameter de- 
pendence. For the boundary condition, the satu- 
ration scale is in the soft (for small va) to semi- 
hard (for large va) regime. We need to go to very 
small x ~ 1CP 6 to finally leave the soft region be- 
hind. Notice however that for most observables 
an average over va is implied, and the contribu- 
tion from the impulse approximation (or lowest 
order in j) is never really small (see e.g. |l|6j ). 
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3. Photon- Jet correlations 



It would now be helpful to find an observable 
which allows to map the p-dependence of the un- 
integrated glue, and to determine the saturation 
scale as directly as possible. We have no space 
here to explain here in detail the formalism of 
nonlinear k± -factorization, which principal fea- 
ture is that cross sections of hard processes in 
a nuclear environment are highly nonlinear func- 
tionals (quadratures) of the collective nuclear un- 
integrated glue [1171618] , The emerging nonlin- 
ear k± -factorization formulas for the hard dijet 




spectra we shown to fall into several universality 
classes depending on color properties of the un- 
derlying pQCD subprocess [21618] . These classes 

differ in the character of incoherent initial state Figure 2. Top:i? p ,4 at x = 0.01, Lower: R cp for 

and final state interactions and coherent distor- v > = 8, v < = 1 for different x. 

tions of two-parton Fock states of the incident 

parton. Of special interest are reactions in which 

both the direct photon and the accompanying 

(balancing, recoiling) quark jet are observed. We 

shall refer to such final states as dijets. The pro- thc prob lcm abelianizes, and all multiparton Se- 

duction of such qi dijets can be viewed, in the matrices can be reduced to the ones of gg-states. 

nucleus rest frame, as an an excitation of the The emerging dijet spectrum is then a /mearfunc- 

57 Fock state of the incident quark, q* -> q~f. tiona i f tne unintegrated glue. 

As with all dijet observables |1|2|6|8], the di- 0n the free nuc i e0 n target, we obtain for the 

jet spectrum is obtained in terms of multipar- parton-level cross section- 

ton S-matrices of four (qjqj)-, three(gg7)- and 
two(gg)-parton states. In distinction to purely 

QCD-processes, however there is a great simplifi- 2(2n) da^iq — » qj) _ ^. ^ ^ . . 

cation due to thc fact that the 7 does not interact dzd 2 pd 2 A. 

through multigluon exchanges. As a consequence, xK(p,p — z A) , (6) 
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with 



^(Pl,P 2 ) = 



Pi 



Pi 



pi 



p\ 



and the splitting function 
1 + (1- 



P iq {z) = 2ei 



zf 



(7) 



(8) 



Here z,p are the photon's light-cone momentum 
fraction and transverse momentum, e 2 = zm 2 . 
The decorrelation momentum A = p + p q = 
measures the deviation of the dijet sustem from 
the back-to-back situation. It is worth to point 
out that the spectrum © is exact over the full 
phasespace of the system. In particular notice 
the collinear pole at p = zA from the final state 
photon emission of the scattered quark. It corre- 
sponds to the "monojet" configuration, where the 
<?7 system recoils against a jet at a large distance 
in rapidity. 

For the nuclear target we obtain: 

{2n) 2 da A {q -> yy) = 
dzd 2 pd 2 Ad 2 b 

(v A ,x,A)+w (v A )6W(A) 

xP 7q (z)K{p,p-zA). 



(9) 



A main result of this work is that the decorre- 
lation momentum distribution maps out the nu- 
clear unintegrated glue. Notice that a potential 
diffractive contribution oc wq{va), which would 
violate the linear /c^-factorization, vanishes on a 
heavy nucleus, where the momentum transfer dis- 
tribution is (5-function-like. For a recent discus- 
sion of 7-particle correlations in the framework 
of the Color-Glass-Condensate model, see [5]. 

Nuclear effects are conveniently characterized 
by the ratio (here we stay at the parton-level 
throughout) 

d 1 a \ da a (j>(fA,x,A.) 

RpA[VA,P, A = = —- . 

T A {b)da N va f{x, A) 

A similar ratio is the central-to-peripheral ratio, 
which involves only nuclear quantities: 

Rc P {v>,v<,P, A) = — — 

v> <f>{v<,x, A) 

Interestingly, both these ratios do not depend on 
the photon's transverse momentum p. We show 



the ratio R p a at x = xa for different opacities 
in the top panel of figlU We observe a shadow- 
ing at small values of A and a Cronin-type peak 
which position reflects the ^-dependent satura- 
tion scale. In the lower panel we show R cp and 
its evolution with x. While it displays the same 
Cronin-peak as R p a for x = xa, the latter is en- 
tirely quenched at small x. 

4. Conclusions 

We presented a nuclear unintegrated glue from 
the solution of a nonlinear small-x evolution 
equation. Photon-jet correlations in the proton 
fragmentation region of pA collisions have been 
suggested as a direct probe of the nuclear unin- 
tegrated glue. They may present an opportunity 
to measure the saturation scale at the LHC. 
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